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Abstract: One-dimensional light harvesting structures with a realistic 
geometry nano-patterned on an opaque metallic film are optimized to render 
high transmission efficiencies at optical and infrared frequencies. Simple 
design rules are developed for the particular case of a slit-groove array with 



C/3 

(— I ■ a given number of grooves that are symmetrically distributed with respect to 

O ^ a central slit. These rules take advantage of the hybridization of Fabry-Perot 

modes in the slit and surface modes of the corrugated metal surface. Same 
design rules apply for optical and infrared frequencies. The parameter space 
of the groove array is also examined with a conjugate gradient optimization 
algorithm that used as a seed the geometries optimized following physical 
intuition. Both uniform and nonuniform groove arrays are considered. 
The largest transmission enhancement, with respect to a uniform array, is 
obtained for a chirped groove profile. Such enhancement is a function of the 
wavelength. It decreases from 39 % in the optical part of the spectrum to 15 
% at the long wavelength infrared. 
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1. Introduction 

Exciting new developments in optoelectronics and chemical sensing are based on a single sub- 
wavelength aperture milled in an opaque metal film, which surface is sculpted at the scale of the 
wavelength HI 121. In systems like the slit-groove array (SGA) sketched in Fig.[T] the surface 
corrugation acts like an antenna to couple the incident light into surface modes that squeeze 
the EM energy into the aperture, leading to high transmission efficiencies at a narrow spectral 
range. Surface plasmon polaritons (SPPs), modified by the metal corrugation, are responsible 
for light harvesting in the optical regime; while in the THz regime, where metals behave as 
perfect electric conductors (PEC) and SPP modes no longer exist, surface modes are induced 
by the geometry fS) . 

Tuning the optical response of single apertures surrounded by periodic grooves are of par- 
ticular interest in the context of spectral and polarizing imaging lH |5] IS El . Transmission res- 
onances at a given wavelength can be controlled by adjusting the size and geometry of the 
SGA. Its optical response has been optimized by varying groove depth H, periodicity 
aperture size [8|, metal thickness 0|, and number of grooves llOI . Such studies are not limited 
to uniform and periodic groove arrays, where all grooves are identical and equidistant, but also 
consider nonuniform and not periodic structures lfTTl[T2l[T3l[T4l[T5] |. 

The correlation between geometrical parameters for an efficient light harvesting process 
(LHP) has been discussed in several of the previous references. Authors of Ref. Q have pointed 
out the role played by groove cavity modes, waveguide (Fabry-Perot) modes in the slit, and the 
in phase groove re-emission (controlled by the period of the groove array). They employed 
the PEC approximation, which has only indicative values at optical frequencies. Janssen et al. 
refined the numerical exploration of the array geometry explaining how to squeeze additional 
light into a slit in a real metal by increasing the number of grooves but reducing the size of each 
groove at the same time ifTOI . 

The aforementioned studies are limited to some relevant geometrical parameters and often 
to a fixed wavelength. Moreover, whereas such studies are mainly focused on optical, THz 
and microwave frequencies, less attention has been paid to the infrared (IR) regime, for which 
one can envisage potential applications in spectroscopic chemical sensing |[T6l[T7l . bio-sensing 
lfT8l[T9ll20l . and light harvesting technologies, see Ref. 1211 for a recent review. In the present 
paper, we go one step further and study a broader area of the parameter space and a larger 
spectral range that extends from optical to IR frequencies. 

Our first goal is to identify simple design rules for optimal light transmittance by the SGA 
shown in Fig. [T] Such rules are developed from the physical mechanisms responsible for 
squeezing light into the central aperture. Our second goal is to further optimize the light har- 
vesting process with help of a conjugate gradient (CG) maximization algorithm 1221 that scru- 
tinizes the parameter space of the groove array. Seeds for the CG algorithm are taken from 
simulations based on physical intuition. We shall consider uniform and periodic groove arrays 
as well as nonuniform and not periodic ones, which further enhances the light transmission 
through the slit. Notice that Ref. lT3l has already explored numerically nonuniform SGAs but 
only in the microwave regime. Moreover, a remarkable difference with previous works is that 
we optimized the LHP for a fixed size of the system. Thus, we perform a full optimization for 
a given number of grooves and a given wavelength. Fixed size optimizations are more fruitful 
for practical applications, like optical detectors, which usually have a given pixel size and its 
fabrication is limited by several geometrical constrains. 

Calculations are done in the framework of the coupled-mode method (CMM) ||2l with surface 
impedance boundary conditions ll23l . Such approximate boundary conditions are applied not 




Fig. 1. (Color). Schematic representation of the slit-groove array on a free-standing gold 
film of thickness hs. The slit of width vvj is surrounded by 2Ng grooves (width Wg, depth 
hg, and period P). The distance from the slit to the first groove, dsg, is in principle different 
from P. 

only at the horizontal air/metal interfaces but also at the vertical walls of slit and grooves. This 
model nicely reproduces experimental results ifTl l24l |25]| . 

Despite the intensive optimization presented here, a more efficient LHP can be achieved 
with slight variations of the geometry of Fig.[T] Such kind of improvements have been largely 
considered in the literature, e.g. the inclusion of Bragg reflectors 1261 . the change of the shape 
of either the aperture 1271 or the grooves l6|, or the combination of single apertures perforated 
on a metal film with dielectric surface gratings |28l; see Ref. Q for a comprehensive review. 
The study of such modified geometries exceeds the scope of the present paper However, we 
hope that the present discussion could motivate further theoretical and experimental studies on 
that direction. 

The paper is organized as follows. The coupled-mode method is briefly described in the next 
section. Section|3]studies a uniform SGA. Its optical response is optimized first by a procedure 
based on physical intuition and next with the CG algorithm. Both wavelength and size of the 
system are varied in the optimization process. Section |4] considers a nonuniform SGA of a 
given size as a function of the wavelength. At the end of the paper our main conclusions are 
summarized. 

2. Theoretical framework: coupled-mode method 

We study the system represented in Fig, [T] It consists in a slit of width Ws perforated in a 
free-standing gold layer of thickness hg. The front metal surface is corrugated with a set of 
2Ng grooves (of width w^ and depth hg) symmetrically distributed with respect to the central 
slit. The groove period is P. The distance from the slit to the first groove, dsg, could be in 
principle different from P. We assume hereafter that the SGA is illuminated by a plane wave 
which electric field is in the plane of the surface and perpendicular to the slit direction. The 
dielectric constant of gold is taken from Ref. ll29l and l30l for optical and infrared frequencies. 



respectively. 

We compute the intensity of the Hght radiated to the farfield and normalize it to the intensity 
of the light incident on the area of the slit. The resulting normalized-to-area transmittance, 
called Tj from now on, accounts for the efficiency of the light harvesting process: 77 is of the 
order of 1 for a single sht, whereas it could become one or two orders of magnitude larger when 
the groove array squeezes additional light to the central slit. 

Maxwell's equations are solved self-consistently using the coupled-mode method ll). It 
is based on a convenient representation of the EM fields. Above and below the metal film 
shown in Fig. [T] the fields are expanded into an infinite set of plane waves with both p- and 
s-polarizations. Inside slit and grooves the most natural basis is a set of planar waveguide 
modes ISTI . Convergence is fast achieved with a small number of such modes |[32l . The parallel 
components of the fields are matched at the metal/dielectric interface using surface impedance 
boundary conditions (SIBCs) ||23l . Although SIBCs neglect the tunneling of EM energy be- 
tween the two metal surfaces, this effect is not relevant for a metal thickness larger than a few 
skin depths. SIBCs are also used at the lateral walls of slit and grooves. After matching the fields 
at the interface we arrive to a linear system of tight binding-like equations, than can be easily 
solved. We direct the interested reader to Ref. Il33l for the expressions of the full multimode 
formalism as well as their derivation, and to Ref. If32l for its application to one-dimensional 
systems. 

3. Uniform slit-groove array 

We consider first a uniform SGA, where all grooves are identical. Section [TTI describes the 
physical mechanisms underneath in the light harvesting process and how to combine such 
mechanisms for having an enhanced transmittance. The optimal geometry obtained in this way 
is used in sections 13.21 and 13.31 as a seed for the conjugate gradient maximization method. In 
section [32]the SGA is optimized as a function of the wavelength for a fixed number of grooves, 
while in section [33] A is keep constant and A'^ is allowed to vary. 

3.1. Optimization following physical intuition 

The normalized-to-area transmittance of a SGA can be enhanced when a Fabry-Perot mode is 
excited inside the slit or when the geometry of the groove array supports groove cavity modes. 
For a given wavelength, the Fabry-Perot mode can be tuned by metal thickness /ij (see IJ) 
and references therein), while the groove cavity mode of a groove array is a function of the 
groove depth and period 1341 . However, the largest transmittance is obtained when the Fabry- 
Perot mode of the slit is located at the same spectral position than the groove cavity mode, 
as pointed out in Ref. fS). We find that this recipe applies not only for optical frequencies but 
also for the IR part of the spectrum. Fig.|2ja) illustrates this behavior for A = 1.35 pLra and 
Ng = 10 grooves. The Fabry-Perot peak for a single slit is centered at A = 1.35 /xm for a film 
thickness of = 360 nm. Despite the intensity of the Fabry-Perot mode changes monotonously 
with the aperture size Ws of a sub-wavelength slit 121 [35)1 . choosing the value of Ws implies a 
compromise between a high normalized transmittance Tj that decreases with Ws and a high total 
transmittance T that increases with Ws- Similar trends are found for a SGA. We therefore keep 
constant the value of w^. We use = 100 nm that gives 77 = 1.8 for a single slit at A = 1.35 

The groove cavity mode coincides with the Fabry-Perot mode forP = 0.93Ajpp and h^, = 120 
nm, leading to the largest 77 = 54 in the black line of Fig. I2a). P is normalized by the SPP 
wavelength (Xspp = 1-33 /xm) for the sake of convenience. P < Xspp due to penetration of the 
field in the grooves that enlarges the optical path of the SPPs. SPPs are emitted in phase when 
Xspp is equal to the sum of P and the extra optical path provided by the grooves. We assume for 
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Fig. 2. (Color), (a) Normalized-to-area transmittance (rj) for a SGA as a function of the 
groove period P for groove deptli increasing from h^=l 10 nm to 130 nm at A = 1.35 flm. 
P is normalized to Xspp =1.33 jlm. (b) T] as a function of the slit-first groove distance dgg 



for several values of P/X^pp 
rj versus the aspect radio (r = 



taken from the black curve of (a) (h„ 



125 nm). (c) 



Wg/hg) for the optimal geometry of (b): P = 0.95Xspp and 



d,g = 0.93P. 



the moment that Wg ~ hg and d^g = P. These constrains are relaxed in what follows. 

Fabry-Perot and groove cavity modes are in principle hybridized in a SGA, but roughly 
speaking we could say that the interaction between slit and groove array can be controlled by 
the distance from the slit to its nearest groove, d^g. Fig. |2tb) renders 77 as a function of dsg 
(normalized by F) for several values of P/Xspp- The groove pitch increases from the value 
P = Q.95Xspp, at which rj has a maximum in Fig.|2](a), to the value P = l.Q3Xspp, at which rj 
has a minimum. Fig. |2](b) shows a double-peak when P < X^pp: one peak at dsg ~ 0.4P and 
another at dsg 0.95P. The double-peak is shifted to smaller values of dsg when P approaches 
Xspp. Only a single broad peak lasts when P > Xspp. 

The highest emittance is obtained when dsg is close to P but always smaller than Xspp; for 
our choose of parameters T] = 62 for dsg = 0.93P and P = 0.95 Xspp, see the black line of 
Fig. lUb). It is worth to stress that a slightly larger P implies a significant reduction of 77. It is 
commonly assumed that the ideal groove pitch for extraordinary optical transmission should 
fulfill SPP constructive interference conditions (P = nXspp, where n is an integer), thought it 
does not occur in practice yj; see Ref. [36J for a discussion in the case of a minimal interactive 
system. For example, Janssen et al. obtained an optimal transmittance for dsg = 0.54P assuming 
P equal to Xspp at A = 800 nm lITOl . However, despite of the disagreement with authors of Ref. 
ifTOl regarding the best choice for P, our calculations confirm their results that the transmittance 
is negligible when dsg = P = Xspp. 

The aspect ratio r = Wg/hg is another important design parameter Fig. |2](c) shows 77 as 
a function of r for the optimal geometry of Fig. |2](b). It has two peaks centered at r « 1.0 
and r « 6.0. The position of the second peak changes as a function of the wavelength (not 
shown). Though 77 is slightly larger for r = 1 .0 than for r = 6.0 in Fig.|2](c), we shall see below 
that the relative contribution of the two peaks changes when several parameters are optimized 
simultaneously. 

We conclude the analysis of Fig. |2] stressing that groove pitch and depth are the most rele- 
vant design parameters, leading to 7j = 54 in Fig. |2|a). Ideal values of r and dsg produces an 
additional enhancement of 15 %. Up to now we have rj =62 for A = 1.35 jim, hs = 360 nm, 
Ws = 100 nm. hg = Wg = 125, P = 0.95 Xspp, and dsg = 0.93P. This geometry is used below 
as a seed in the CG optimization. Similar calculations have been done for other wavelengths at 
optical and IR frequencies (not shown). 

3.2. CG optimization as a function of the wavelength 

We search now for the optimal groove profile scanning the parameter space with help of the 
conjugate gradient maximization algorithm 1221 . that was implemented for that purpose within 
our CMM formalism. Optimal geometries as a function of A for a uniform groove array with 
Ng = \Q grooves are reported in Table [T] 

The metal thickness (hs) defining the position of the Fabry-Perot mode for a given wavelength 
and the slit width (w,) that defines its intensity are keep constant along the simulation process. 
For the sake of comparison on an equal footing, we scale Ws at a given A with respect to the 
value Ws = 100 nm chosen above for A = 1.35 /im. Input values of hs and Ws are included in 
TablelU 

We observe in Table [T] that the normalized transmittance for the optimal system increases 
from 7j = 35 at A = 0.65 nm to 7j = 74 at A = 10.6 /zm, due to the reduction of the metal 
absorption for increasing A. The absorption is also given in Table [T] It decreases from 59 % at 
A = 0.65 /im to 22 % A = 10.6 /im. This quantity is defined as 1-R-T, where R and T are the 
total reflectance and total transmittance, respectively. 

The most remarkable feature is that, despite the reduction of the absorption, 7j does not 
increase with A in the mid-IR but states almost constant around the value 7] = 75. That indicates 
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Table 1 . Optimal geometries obtained by a CG optimization of a uniform SGA with Ng — \Q 
grooves for different values of the wavelength. Last column reports calculations for a PEC 
at A = 10.6 lim. Last row gives the values of r\ for a SGA scaled by X with respect to the 
SGA made of gold and optimized at A = 10.6 /im. 



that radiative losses dominate over dissipative losses. However, Table [T] shows that dissipative 
losses are not negligible even in the IR regime. In order to quantify its effect we optimized a 
SGA in a PEC at the longest wavelength considered in this paper, X = 10.6 jim. Results are 
given in the last column of Table [U We find that, even at the long wavelength part of the IR, T] 
for a PEC is 42 % larger than for gold. A relevant feature of a PEC is that the dimensions of 
the system can be scaled with X, having the resulting structure the same transmittance than the 
original one. We have verified if that simple rule could be applied to the IR. Our calculations 
show that it is a nice approximation for the long wavelength part of the IR. In the last row of 
Table [U we report the values of 77 for the SGAs scaled by A with respect to the SGA made of 
gold and optimized at A = 10.6 jim, i.e. we multiply the geometrical parameters of the before- 
last column of Table [T]by the scaling factor 8^, which is equal to 8x = 8.0/10.6 = 0.75 for 
A = 8.0 /im, 5^ ~ 0.57 for A = 6.0 /xm, and so forth. A good agreement between the optimal 
SGA and the scaled one is obtained for A > 6.0 /im. 

We find however that the ideal P « Q.95Xspp and d^g ~ 0.92P are independent of the wave- 
length for A > 0.65 /im. These values are very close to those obtained for a PEC. This behavior 
is related to the properties of groove cavity modes ll32l . The dispersion relation of surface 
modes is close to the folded dispersion relation of the flat surface, which is also close to the 
folded light line (notice that Ajpp w A in Table [T]i. At the center of the Brillouin zone, gaps are 
open due to the finite size of the grooves. Gaps strongly depend on groove geometry and wave- 
length, leading to an ideal aspect radio that decreases from 6.6 to 2.6 for increasing A, as well 
as an optimal groove depth that oscillates around a quarter of the metal thickness. Geometrical 
changes are larger for a PEC due to the different scattering cross-sections of real and perfect 
electric conductors. 

The tolerance of the optimal system to small changes of its geometry is estimated by the full 
width at half maximum (FWHM) of 77 as a function of the corresponding geometrical parame- 
ter. The FWHM for the optimal geometry given in Table[T]is reported in Table|2l In particular. 
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Table 2. FWHM of r] as a function of each parameter reported in[T] In particular, the 
FWHM for ri versus X is the bandwidth. Last column reports calculations for a PEC at 
1 = 10.6 lim. 

the FWHM for rj versus A is the bandwidth. We find that the tolerance on the period is smaller 
than for dsg- The FWHM is also larger for Wg than for lig. The reason for this behavior is that 
period and groove depth give rise to the largest enhancement of Tj when groove cavity modes 
are excited. So, peaks are narrower when the system is more resonant, as one should expect. 
Interestingly, rj versus Wg has a broad peak (not shown) instead of the two peaks depicted in 
Fig. He). 

5.5. CG optimization as a function of the system size 

In this section we optimized a uniform SGA as a function of the number of grooves, A''^. A 
CG maximization is performed for a given value of Ng and A using as a seed the geometries 
reported in Table[T] Fig. [3]shows our results for Ng increasing from 2 to 42 and A = 1 .35 /im. 

We observe in Fig. [Sja) that the normalized-to-area transmittance grows monotonically up 
to Tj =216 for Ng —42. As pointed out in Ref. [[TOl . in order to achieve such enhancement the 
distance between grooves raises as a function of Ng and the groove depth decreases, see Figs. 
[3](b) and (c). In particular, we find that P approaches Xspp for increasing A^^. In contrast, dsg 
reaches its maximum value for only 14 grooves. It means that the central slit no longer interacts 
with the outer grooves when A^^ is larger than this value. On the other hand, hg decreases by 
100 nm as a function of A^^ until it reaches the value hg = 72 nm for Ng = 22. Then, hg remains 
practically constant. We add to the conclusions of Ref. iflOl that, in order to perform the optimal 
LHP when more grooves are added to the SGA, grooves should be not only shallower but also 
broader, see Fig.[3](d), where the aspect ratio raises until r = 11.6. Notice that only Wg still 
grows for A'^ > 20 whereas P, dsg, and hg are practically constant. We also see in Fig [3a) 
that the FWHM decreases with A'^. Thus, the price for the enhancement of 77 is having a more 
resonant SGA. 

Features found in Fig.|3]can be explained as follows. Given that the light harvesting process is 
mainly dominated by radiation losses, the SGA increases its total size to collect additional light, 
thus increasing P and d^i,. For the same reason, broader grooves enhance the scattering cross- 
section between modes inside the grooves and the EM radiation |34J. On the other hand, as the 
radiation pattern is driven by surface modes running along the metal/air interface, the grooves 
tend to be shallower to reduce the shadows that apertures closer to the central slit produce on 
apertures on the edge of the SGA. Shallower grooves also implies a smaller scattering cross- 
section, but this effect is compensated with broader grooves and a larger system size. 
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Fig. 3. (Color). CG optimization of a SGA for A = 1.35 and A'^ increasing from 2 to 
42. (a) Normalized-to-area transmittance, rj, and FWHM. (b) Period, P, and slit-nearest 
groove distance, dsg- (c) Groove depth, hg. (d) Aspect radio, r = Wg/hg, and groove width. 



4. Nonuniform slit-groove arrays 

We can also vary the geometry and pitch of all the grooves, in such a way that the optimized 
groove array becomes aperiodic and the groove geometry change with the groove position. In 
such a case, physical intuition can not longer be used and the optimization process is based 
solely on the CG algorithm. The optimization routine is divided into the following steps: 

i) Optimization of a uniform SGA: already discussed in the previous sections. 

ii) One-parameter optimization: A given parameter is let to vary for each groove, so the groove 
array is no more uniform. Independent simulations are run for the set of groove depths, {hg}, 
groove widths, {wg}, and groove positions, {xg}. 

Hi) Two-parameter optimization: A pair of parameters is let to vary for each groove. Indepen- 
dent simulations are run for the sets: {lig,Xg}, {lig,Wg}, and |wi^,,X|^,}. 

iv) Full groove geometry, {hg,Wg,Xgj: All geometrical parameters in each groove are let to 
vary. 

Optimal geometries obtained in one step are used as a seed for next step. Several trials with 
different seeds have been done for checking the consistency of the numerical results. We find 
that the variation of each parameter is within the tolerances reported in Table |2] Spectra for 
SGAs optimized at several Xs are represented in Fig.|4] Typical profiles are depicted in Fig.|5] 

We find that the enhancement provided by a nonuniform SGA (optimizing the full groove 
geometry) with respect to a uniform one is 39% at A = 1 .35 jim, 17 % at A = 4.0 /im, and 15% 
at A = 10.6 jim, see Fig.|4] The enhancement decreases with the wavelength. 

The variation of the groove depth provides optimal geometries among one-parameter opti- 
mizations, with an enhancement of 11 % for A = 1.35 /im, 10 % for A = 4.0 /im, and 8% 
for A = 10.6 jim. Fig. [Sja) shows the groove profiles obtained at A = 1.35 fim after the 
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Fig. 4. (Color). Spectra of the SGA optimized with the CG algorithm for the target wave- 
lengths (a) A = 1.35 flm, (b) A = 4.0 jlm, and (c) A = 10.6 ^m. Full {hg,Wg,Xg} and 
one-parameter optimizations ({/ii;}, {wg}, and {x^}) are compared with the optimization 
for a uniform SGA. 



optimization of a uniform SGA (red), the one-parameter optimization j/ig} (black), and the 
full-parameter optimization {hg,Wg,Xg^ (blue). The optimal groove-depth profile for {hg} is 
a chirped array with groove that deepens from 0.24/1^ to 0.58/1^. A similar variation between 
O.2I/I5 to 0.64/2s is observed for the full parameter optimization. We also find that grooves oc- 
cupy a smaller fraction of the metal thickness for longer wavelengths, see Fig.|5jd). In all these 
cases, distant grooves need to be deeper than those close to the central aperture in order to 
increases their scattering cross section and squeeze additional light into the aperture. 

The optimal width profile is also a chirped array with incremental groove width, see Fig.|5lb). 
As for the depth profile, distant grooves are broader for they should increase their scattering 
cross section in order to further enhance the transmittance through the central aperture. We also 
observe a difference between {wg^ and {hg,Wg,Xg^ optimizations larger than for the groove 
profile. The optimal geometry is consistent with previous calculations that found an efficiency 
of 50% in plasmon generation using a grating with a gradient in the groove width fTZl . 

We do not find large deviation from the uniform SGA when the groove-position profile is 
optimized with the CG algorithm, see Fig.|5|c). The relative distance between grooves is prac- 
tically the same than for a uniform SGA, except for the farthest groove of the array in the 
full-parameter optimization. The optimization of the groove-position profile does not produce 
an appreciable enhancement of the normaUzed transmittance with respect to the uniform SGA, 
see Fig.E] 

5. Conclusions 

In search for high efficiencies in the light harvesting process of a one-dimensional slit-groove 
array, we have considered the interplay between Fabry-Perot modes of a single slit and groove 
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Fig. 5. (Color). Profiles obtained for A = 1.35 flm in the CG optimizations of the uniform 
and nonuniform SGAs. (a) Profiles of groove depths, (b) Profiles of groove widths, (c) 
Profile of groove positions, (d) Groove profiles as a function of A. 



cavity modes of the groove array. We have developed the following simple design rules for a 
uniform SGA, that are valid for optical and IR frequencies. 

• Fabry-Perot and groove cavity modes should be at the same spectral position. For sub- 
wavelength apertures at a given X, the position and intensity of the Fabry-Perot is con- 
trolled by metal thickness and aperture size, respectively. The spectral position of the 
groove cavity modes is mainly determined by groove depth and pitch. 

• The normalized-to-area transmittance grows monotonically with the number of grooves, 
given that grooves move apart and become broader and shallower as A'^ increases. In 
particular; 

- The ideal periodicity (?) is slightly smaller than the SPP wavelength (Xspp). P tends 
to Xspp (from below) when more grooves are added to the SGA. 

- Slit-nearest groove distance should be about Q.92P (but only for the optimal P). 

- The ideal groove depth is smaller than half of the metal thickness. This quantity 
decreases with the number of grooves (A'^;.)- 

- The optimal aspect ratio decreases with X but raises with Ng. 

We also account the following trends for the optimal system: 

• The transmission efficiency of a uniform SGA in the infrared is practically independent 
of A. 

• In contrast, the enhancement provided by a nonuniform SGA decreases with A . A chirped 
groove array enhances the transmittance between 15% and 39% for decreasing A. 



• The bandwidth decreases with the size of the system. 

Summing it up, we hope that our findings could motivate further theoretical and experimental 
studies of light harvesting structures. 
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